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Abstract 

I We suppose and develop a simple quantitative model of polymer film conductivity. This model 

can be seen as a further development of the ideas of Vlasov, Apresyan et al pQ. The main point of 
the model is that conducting islands exist, and the charge transfer between the islands is carried out 
by mobile segments of polymer molecules. This model quantitatively describes the presence of two 
O I states of conductivity, and the current stabilization phenomena, and it predicts the temperature 

dependence of conductivity, and the dependence of conductivity on the thickness of the polymer 
film. The pressure-driven transition to high-conductivity is described only qualitatively in this 
<^> model. 
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1 Introduction 

1.1 Motivation 

In recent years interest in the conductivity of polymer films has greatly increased [2], [3] , and several 
Nobel prizes in physics and chemistry were awarded in recognition of people's success in polymer film 
conductivity. There is a plethora of possible applications of polymer electronics in different fields: in 
science, in technical devices, in computers, etc. But to use polymer electronics we should understand 
polymer conductivity better, or at least create a model which answers the most simple questions: 
"Why and how do the polymers conduct current? What are the two states of conductivity? How can 
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we switch between them and why? Why does the stabilization of current occur?" and so on; there are 
a lot of questions. This paper pretends to be a toy model of polymer conductivity, which collects and 
explains all notable facts known about the conductivity in plasticized PVC (polyvynilchloride) : the 
two states of conductivity and quasi-spontaneous transitions between them the current stabiliza- 
tion in the state of low conductivity [5], the pressure-caused transition to a high conductivity state, 
and the relaxation of current upon the constant voltage. In addition, the model gives the quantitative 
answers to all the above-mentioned questions, excepting the pressure-induced transition to high con- 
ductivity. The pressure-induced transition is explained only qualitatively in this model. This model 
was developed by using the facts about the conductivity of PVC films, but it is supposed to catch 
some general patterns of a great variety of different polymers. 

1.2 Outline of the model 

Our model is a development of the ideas of D. V. Vlasov, L. A. Apresyan et al p]. The main ideas 
are the following. The structure of polymers is known to be very complicated and involute, sometimes 
exhibiting fractal behaviour. So, the density in polymers is not a constant, but it varies from place to 
place and depends on the characteristic size on which we calculate the density. In the regions of high 
density the atoms of different molecules can approach each other so closely that the wave functions of 
electrons begin to smear between two or more different atoms. In such regions the delocalisation of 
electrons can occur and thus the electrons can move within such an area almost freely. This area of 
a polymer may be similar by its conducting properties to semiconductor. Further we call such areas 
"conducting islands". The only remaining question is - how is the charge transferred between the 
different conducting islands? In this model the moving segments of polymers molecules are believed 
to do it, see Fig{T] Without plasticizers the PVC molecules are like rigid sticks, but plasticizers allow 
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Figure 1: The scheme of charge transfer between islands. 1,2 - conducting islands, 3 - moving molecule 
segment, 4 - the movement of the segment 

them to bend. Some of such segments are within reach of two islands, as shown in the picture FigjTJ 
The sequence of conducting islands, connected with such junctions, we call the conducting channel. 
The segment can have a defect, which can hold some charge, most probably several unit charges. This 
defect can carry both signs of charge, but a preferred sign may exist. If so, there are two modes of 
conductivity. In the first mode, the segment moves in one direction with a charge, but in the opposite 
direction without a charge, only due to temperature fluctuations. Furthermore, the second mode is 
where the segment moves with a charge in both directions. When the segment moves with a charge, it 
moves upon action of the electric field, and therefore it moves faster, given that field is more than some 
critical field. These two modes of conductivity are the two conductivity states, and if the movement of 
a segment upon the field is faster than the movement upon the temperature fluctuations, the current 
in low-conductivity state is roughly constant with applied voltage, as it was observed in experiments. 
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The transitions occur upon the high voltages, and upon suddenly supplied voltages. 



1.2.1 Pressure-induced phenomena 

The pressure can have the same consequences as the applied field: due to the pressure, the neighbouring 
conducting islands become closer, and the field between them increases. In addition, the pressure 
causes the charge division even without the external field. This is because the reaction force which 
appears when the pressure is supplied is an electromagnetic force (not gravitational, weak or strong), 
hence, the charges should exist to produce electric force. This charge division in polymer materials 
can be easily observed: your plastic comb attract small pieces of paper after you have used it. The 
appearance of the macroscopic force tells us that huge charges play a role in this phenomena - this 
experiment can be used to reinforce the argument that pressure causes great charge division. On the 
other hand, the applied field causes the charge division, and at different sides of each junction are the 
charges of different signs. These charges attract each other, creating the pressure in the junction. 

2 The description of the model 

2.1 The dependence on the thickness of the polymer 

In our model, it is supposed that the current goes by some channels which are nothing but the system 
of islands, connected with a conducting junction. Indeed, not any junction can conduct current. To 
conduct current, junction must have an appropriate mobile segment, and this mobile segment should 
be appropriately oriented. The existence of mobile segment itself require the existence molecule of 
plasticizer in junction, because it is the plasticizer molecule which enables the PVC molecule to be 
not rigid. From this we can infer, that the more the concentration of plasticizer is, the more the 
conductivity is. More detailed consideration is postponed to Sect |2.5| Additionally, these molecules 
should be properly oriented. From all these considerations we can infer, that in our model not each 
junction can conduct electricity, but electricity flows by some sequences of islands and conducting 
junctions, which we call conducting channels. But the idea of channels itself can tell us a lot of 
information about the dependence of conductivity on thickness. Let's denote the average number of 
conducting channels per unit area in the polymer as f(D), where D is a thickness of polymer film. 
Now we can consider two adjacent films, or, better, one film divided by plane, parallel to the sides of 
the film. The two layer of the polymer have thickness of D\ and D2 correspondingly, D\ + D2 = D. 
Denote the size of conducting island d, and consider the piece of the polymer of area S. We assume 
here that the width of channel equals to the size of conducting island. There are Sf(di) channels from 
the upper layer, and Sf(d2) channels of bottom layer, but the number of channels which go through 
all composite film is lower, because not all channels of upper layer meet channels of bottom layer. Now 
we estimate this number, assuming that 1) the channels are rare 2) they are independent. Consider 
some particular channel of upper layer. To meet with it, the center of bottom layer channel must fall 
into the circle of radius d around the center of the channel under consideration. This region has an 
area of ird 2 . Thus, the number of bottom-layer channels falling in this region is nd 2 /(Z^)- So, one 
channel of upper layer increases the number of channels through entire film by ird 2 /(D2), and there 
are Sf(D\) upper layer channels. So, the number of channels in entire film equals: 



The only arbitrary constant here is do, the coefficient in front of exponent is fixed by the equation. 
Of course, this consideration is rather robust, it neglects that the channels are of different sizes, that 



Sf(D) = ird 2 f(Di)Sf(D 2 ) 



(1) 
(2) 




This functional equation on f(D) has a following solution: 




(3) 
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they are not independently distributed over area, the layers may not be independent, the edge effects 
etc. 

2.2 The one-channel current 

In our model, there are conducting islands and the charge is transferred between different islands by 
moving segments. Let d be the size of conducting islands and cxo be their conductivity. Of course, 
both the size, conductivity and other parameters are different for different islands, d and <7o represent 
the average values of these parameters, and this simple model may catch the general patterns of the 
conductivity. Denote the length of the mobile segment as I, and an angle it turns by as a (see Fig§). 
The linear density of mobile segment is pi - but it is not a free parameter of the model: it is determined 




Figure 2: The parameters I and a are depicted on this figure. 



by the structure of polymer molecule. Denote the dimensionless parameter 7 = d/{2 sin(a/2)Z) -this 
is a ratio of the size of islands to the gap between them. The current through one junction equals: 

hjunc = -7-: — , * , - low conductivity (4) 

L\It + l±Ie 

2e e 

hjunc = A , = -r— - high conductivity (5) 

where Atx is time of moving of mobile segment from one island to another upon temperature fluctu- 
ations, AtE is a time of moving of mobile segment from one island to another upon the field acting 
on the carried charge. The current flowing through one island equals: 

hisiand = E d aod 2 (6) 

where E d is a field inside the conducting islands, and E[ is a field between the islands, d 2 in order of 
magnitude equals to the cross-section of the island. In the stationary case, that is when the charge 
density does not depend on time everywhere, the hjunc = hisiand- Then: 

hisiand = hjunc (7) 
Ed&od 2 = hjunc (8) 
E m ac(l + 7) = lE d + Ei (9) 

Here E mac is a macroscopic average field. Once we have found the one-junction current, we can obtain 
the overall current density: 

j = hisiand f(D) (10) 
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Now we are going to discuss the dependence of At? and At e on the parameters of model and applied 
field. Moving segment consists of thousands and possibly millions of atoms, so we can describe it 
classically. Then, its rotation has the mean kinetic energy equal to \kT. So (here I is the momentum 
of innertia of the segment, I = hpil 3 )- 



1 Ico 2 = l kT (11) 
2 2 v 1 



This formula does not take into account the fact that the motion of the segment is rather chaotic 
when moving upon the action of temperature fluctuations. At least, it can undergo the collisions with 
other molecules. If it has collided iV times with other molecules, its mean velocity is decreased by 
There are many ways to estimate N. By order of magnitude, N may be equal to: 

iV=(^-V, (13) 

\ Icon / 

where l con is length of one connection between atoms in molecule. This estimate can be derived as 
follows. During the motion the end of segment will get over the path of al. On this path, it will each 
l con meet with other molecule - so the total number of collisions can be estimated as The motion 
is two-dimensional, so we can take this into account by squaring the expression. So, the final formula 
for Atx is: 

In order to obtain the expression for AtE we use the equation for angular momentum (/ here means 
the moment of inertia of the segment, I = ^pil 3 ): 

^ , 1 + cos(a/2) /q\ , . 

ICj = Eid K 1 ' = eEil cos y-j (15) 

iwAti = a (16) 



At, = (17) 
y 3eEicos (f j 

Since the initial velocity is supposed to be zero, the time of movement is in inverse proportion to the 
square root of acceleration, i. e. to the square root electric field between islands, E[. Note however, 
that this formula is true when the velocity obtained due to temperature fluctuations is negligible 
comparing with the velocity obtained due to the electric field, that is, At? ^> Ai#. Otherwise, 
temperature fluctuations play sufficient role even when moving upon the field, and this can affect the 
resulting time of movement. 

2.2.1 Voltage-current curve in high conductivity 

In this case, the current through one junction is determined by the formula it equals: 

11 = 2A7^ = -KT E (18) 

We write equations (8|9) in this case: 

Emac(l + 7) = lE d + Ei (19) 
a E d d 2 = (20) 

L\lE 
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3e 3 cosf — ) 

We introduce here the new derived parameter, Eq = 2aa i d^pi 2 ' Using the Eq, one can obtain such an 
equation and the solutions of this equation: 



/ 2. / _ ( X + l)Emac . _ / El 

e mac X -\- X , S m ac o TP ' X \ ,2 it* V^-V 

7 Eq y 7^-&o 

e mac = ' E maccr itE = ~ ■ Eq (22) 

^maccritE 1 + 7 

^ = ~ 1±V/1 9 + 4eU: = £(>A + 4eU:-l) - { ^ ''H J ! (23) 

* * IV mac i ^mac ^ 1 

Let's consider first the case when -Emac <C E maccr itE' 



x ~ e mac (24) 

V- 1 - T // -^77 

7 2 



^= V V T C (25) 



j2 l + 7 



Jl = aod E mac (26) 

7 

7 + 1 / £)\ 

j = o" exp - — £" mac (27) 

77T V d 0/ 



o-' = oq^^- exp ( -— ) (28) 



In this limit case, the current is proportional to the applied field. The overall conductivity in this case 

r —exp( D 
77T V d , 

In other limit case, E mac > E maccrtE : 

x = \/ e 'mac (29) 

E l = (j+l)E mac (30) 

h = (Torfy/il + ^EoEmac (31) 

j = ^-^{l + ~f)E E mac exp (~^) (32) 

In this case, the current is proportional to the square root of the applied voltage. This has very simple 
physical interpretation: the current in conducting islands is proportional to applied field, and the 
current through junctions is proportional to the square root of applied field. These "resistors" are 
connected consequently, so the overall current is determined by the lowest current. Upon low values 
of field, the linear function is smaller, and upon the large values of field, the square root of field is 
smaller. 



2.2.2 Voltage-current curve in low conductivity 



Low conductivity state is a state where the mobile segment moves in one direction due to the temper- 
ature fluctuations, and in other direction it moves upon the force exerted by E\. So, we substitute in 
the expression Q) the formulas ( fl2fl7| ): 



aod 2 



At T + 



2a P il 2 
3e£ { cos(f) 



(33) 
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We denote 



= (34) 

3e 2 cos(f) V ' 

e - ^ (35) 



X= i^k (36) 



and obtain the equation: 



e! 

^"mac 



£X + X + x(l £€ m ac) e mac ~ 0> e mac — (1 + 7) j2jTJ^ (3?) 

Let's solve this equation in two limits. The first limit we consider is £e' mac S> 1. Then the equation 
becomes quadratic: 

ee' mac = 1 + ex 2 (38) 

- = V / <^ ( 39 ) 
Ei = E mac (l + 7 ) (40) 

^ = ^iee' mac J^ (41) 

Usually, 7 > 1. So, when > 1, AtT 3> Aig and the current does not depend on the applied 

voltage - and this indeed is a stabilization of current, which were discovered in the work [5]. In this 
limiting case, 

/ QZ-T 

(42) 

3T = hjuncf(D) = — ^ \l ~vf~n ex P ( ~~T ) ( 43 ) 




The nearly constant curr ent in this case is proportional to V? 1 - this can be verified by experiments 



From the formulas (5|4) one can derive, that = 3 ^ ah , where the j^g/i is a current in high- 
conductivity state, and ji ow is the current in low-conductivity state - so in this limit, the currents in 
the high and low conductivity states have different orders of magnitude. The opposite limit, £e! mac <C 1, 
corresponds to the case when AtT < Aff; - and, as we have already mentioned, this case requires 
more accurate treatment. In this case the time of moving upon the field is not equal to Atg, defined 



by (17) - in this process we should take into account the temperature fluctuations when moving upon 
the field. Let's calculate the limiting field E mac , above which AtE <C Atx- To estimate this, we can 
formally apply the formula Q in the opposite limit case, ee' mac <C 1. In this limit, AtE 3> Aiy and 
the formula Q reduces to the formula So, we obtain that 

1 + T 

h = a d 2 E mac (44) 

7 

We obtain the critical field if we will write I\ = lyr- The critical field value equals: 



7 e / 3kT 
E maccritT - ^ ^ W ^3 (45) 
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Figure 3: The scheme of conducting islands. The left picture is before the transition of middle junction 
in high-conductivity state, and the right is after such a transition. After transition the field in the 
same junction falls in magnitude, but in neighbouring junctions is increases. 



2.3 Transitions between the conductivity states 

The low conductivity state is the state where the segment moves in one direction upon the field, and 
in other direction upon the temperature fluctuations. The high conductivity state is the state where 
the segment moves in both directions with carried charge, in both directions upon the field. How and 
why the transitions between these two states occur? We can name the two reasons for that. The first 
reason was already mentioned in Introduction. When the external field is applied, the charges inside 
the conducting islands are divided, and they are found at the edges of conducting islands, creating 
the field in junction. The charges of different sign are found at the opposite sides of junction - and 
they attract each other, creating pressure in the junction. Due to this pressure, the charge trapping of 
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both signs becomes easier. This also explains, why the effect of applied voltage is partially equivalent 
to the effect of external pressure - pressure can also drive the polymer to high-conductivity state. The 
second reason is that upon big field, the segment with charge acquires big velocity moving through 
the junction. With this big velocity, it will collide strongly with the conducting island. So the more 
the field is, the stronger the interaction between the segment and the island is and the higher the 
possibility to catch the charge is. Catching charge in both mechanisms is a random process which has 
sufficiently non-zero probability in wide range of applied fields, and in both mechanisms the proba- 
bility of catching increases with the field. We can consider some critical field - if the field between 
islands is greater than the critical, the segment will inevitably catch the charge. So, when the field 
in junctions is less than critical, both states can exist, but if the field is above critical field, only 
high-conductivity state can exist. At the same time, the factor exist which makes both states of low 
conductivity and of high conductivity more stable. Consider the polymer which is in low-conductivity 
state. This however does not mean that the field in all junctions is below critical field. Assume the 
field between some two conducting islands is larger than the critical field, then this junction makes a 
transition in high- conductance state. This transition changes the conductivity of junction by orders of 
magnitude, and this causes the rearrangement of field and charge patterns, this is shown at the Figj3j 
One can see, that after transition to high-conductivity state in one junction, the field in this junction 
is decreased, and in neighbouring junctions is increased. This may have different consequences. The 
neighbouring junctions may also make a transition to high-conductivity state, and they will increase 
the field in the next junctions, and this causes them to go in high-conductivity state - eventually the 
entire polymer will go in high-conductivity state. But this excitation may be dumped and may die 
without macroscopic consequences. This may occur because after the junction go to high-conductivity 
state, the field in this junction is decreased, and this junction will tend to go back, in low-conductivity 
state. Due to the same reasons, if all junctions are in high-conductivity state, the transition of few of 
them in low conductivity is unstable, because their conductivity falls, and the field in these junctions 
increases, causing them to go back in high-conductivity. So, we have shown, that both low-conductivity 
and high-conductivity states of entire polymer are locally stable, i. e. are stable to transition of few 
junctions in other state. So, in some range of supplied voltage the both states of conductivity can exist 
in polymer, and the transitions between them are rather stochastic. The lower bound of this range is 
the voltage, upon which in low-conductivity state the field between some islands is more than critical 
field. The higher bound of this range is the voltage, upon which in low-conductivity state the field 
between alj^] islands is greater that critical. If the voltage is greater that the just-named upper bound, 
then the low-conductivity state is globally unstable: if it would exist some time, it would instantly go 
to high-conductivity state, because the field between all islands is greater that the critical field. The 
region of voltages where both states can exist is extended by VI characteristic of voltage supplier. In 
the great majority of cases, the supplied voltage is decreased as the current is increased. This enlarges 
the instability: if the high-conductivity state suddenly appears, the current is increased, the voltage 
falls down, and this encourages the polymer to go back in low-conductivity. 

2.4 Non-stationary case 

All the above analysis has been made in the stationary case. The consideration of non-stationary 
phenomena can extend the range of voltages where both states can exist even more. Consider, for 
example, the suddenly applied voltage. Since the conductivity of conducting islands is higher, at the 
first moment the charges are divided inside them. In any non-stationary case (including this) the 
equation ([7]) breaks down, in this particular case at the first moment Iusiand Iljunc- This creates 
big charges at the sides of islands, and the E\ value can be much more that E\ value calculated using 
stationary formulas. This E\ can become greater that Ei cr i t , so the transition to high-conductivity 
state can occur with very low voltages, but supplied suddenly. The sudden charge dividing in the 
islands just adjacent to the electrodes creates an addition to the current, which looks like a falling 

2 Or, at least, between the great majority 
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exponent and thus it disappears after some time - and very similar relaxation processes were observed. 
The characteristic time of this process is nothing else that the RC time of a single island, whereas 
the value of the current itself depends on the concentration of islands. iiC-time of conducting island 
equals r = RC = ~, where Eq is the electricity constant. The density of current in single island equals 
ji = co-E-macexp (— ^) , so the macroscopic mean density of current is 

jraac = nd 3 (T E mac exp \~~) ( 46 ) 

where E mac is the macroscopic field and n is a concentration (unit per cubic meter) of conducting 
islands. All the islands adjacent to the metal electrode contribute to the falling addition of current, 
not only those who are beginning of conducting channels. In the islands, which are not the beginnings 
of conducting channels, the established electric field equals to E mac , so it appears in this formula. 



2.5 The dependence on the plasticizer concentration 

The concrete form of this dependence can be modelled by some computer experiments. These com- 
puter experiments can give us hints about the structure of the polymer. Assume that the concentration 
and arrangement of conducting islands do not depend on the concentration of plasticizer. Then the 
number of conducting channels is determined only by number of "percolations", that is, the number 
of paths going through conducting islands, connected with each other by appropriate junctions (that 
is, the junction with the plasticizer molecule and with appropriate oriented PVC molecule). The 
concentration of plasticizer determines the probability of two neighbouring conducting islands to be 
connected by appropriate junction. The field along the percolation is approximately the voltage di- 
vided by the length of percolation. Some answers about the dependence of the number of percolations 
on the probability of existence of appropriate junction between two neighbouring islands are given by 
the percolation theory, but in this case what do we need is the number of percolations. In different 
conducting states, we may need a number of percolations weighted with different functions of perco- 
lation length. For example, in low-conductivity state, the one-channel current does not depend on the 
applied field if it is above E maccr itT, so we should take into account all paths with the length less than 
some fixed length. In high-conductivity, we should sum over all percolations with weight of inverse 
length of percolation, etc. So, strictly speaking, the dependencies of conductivity on thickness in 
low-conductivity and in high-conductivity states are different. These dependencies can be calculated 
in computer models which address percolations. 



3 Conclusion 

We propose a quantitative model of polymer conductivity which describes a majority of observed prop- 
erties of conductivity in the PVC films: two states of conductivity, spontaneous transitions between 
them, the stabilization of current in low conductivity state, the relaxation of current upon the constant 
voltage, the pressure-induced transition to the high-conductivity state. The numerical expressions for 
most of these phenomena are obtained, except for the pressure. Pressure-induced phenomena are 
described only qualitatively in this model. The two mechanisms of charge transfer play a role in 
our model: the charge transfer inside the areas where the delocalisation of electrons occurs, the so- 
called conducting islands, and the charge transfer by mobile segments of PVC molecules, which were 
made non-rigid by plasticizer molecules. The conductivity inside the conducting islands is described 
phenomenologically by introducing their conductivity <jq. The charge transfer between the islands is 
described as the charge transfer by the mobile segment. In the low conductivity state, the mobile 
segment catches the charge of one sign and moves upon the field, and it moves in another direction 
upon the temperature fluctuations. In the high conductivity state, the segment moves upon the field 
in both directions, because it carries the charges of both signs. All quantitative formulas are found in 
the Sect j2j but here we collect all the most important dependencies of experiment-observed quantities. 
The stabilized current has the following dependencies: 
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i. j T oc Vt 



2. j T ocexp(-g) 

3. 2t increases with the plasticizer concentration 

There are some directions for the future research. At first, all these formulas are to be compared 
with experimental data. Our model has five parameters: I, d, ao, do, a, but predicts a lot of quantities. 
Secondly, the non-stationary phenomena require more accurate quantitative treatment. For example, 
the theory should predict the period and spectral composition of relaxational oscillations [6], which 
are generated by the polymer film. The dependence of frequency and spectrum of these oscillations 
can provide a lot of data to be fitted with these five parameters. Another direction of possible research 
is the extending of considerations of Sect |2.5| the number of percolations, the weighted number of 
percolation, the dependence of percolation number on the probability of single junction to conduct - 
all these quantities can be obtained from the numerical experiments. 

3.1 The brief summary of formulas 




Figure 4: A sketch of dependence of j on E mac , which is proportional to the supplied voltage. The red 
curve is for high-conductivity, the black line is for low conductivity. The region where E mac < E maccr uT 
and hence AtE > Atx is marked by question sign, because this case requires separate treatment. 



The formula for the stabilized current in low conductivity, Jt, is (43): 



3kT 



JT 



card 2 y Npil 3 



cxp 



D 

do 



(47) 



But if the E mac <C E maccr i t T, this model does not give answer in both states of conductivity, because 
in this case we should somehow make an estimate of the time of movement of mobile segment upon 
both temperature fluctuations and small electric force acting on it. The formulas for current in high 

J I I I 2 

conductivity are (27), (32) depending on the relation between E mac and E maccrit E = yt—Eq: 



Je 



7+1 

oo-^rexp 



D 



do I Emac - 



E„ 



-'mac ^ E maccr itE 
^ a/(1 + ^EoEmac exp f — lr J , E mac » E maccr i t E 



(48) 



All these relations are schematically shown at the FigJIJ 
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